Combined Influence of Off-diagonal System Tensors and Potential 
Valley Returning of the Optimal Path 

Chun- Yang Wang£ 
Department of Physics, Shandong Provincial Key Laboratory 
of Laser Polarization and Information Technology, 
Qufu Normal University, Qufu 273165, China 

Abstract 

The two-dimensional barrier passage is studied in the framework of Langevin statistical reactive 
dynamics. The optimal incident angle for a particle diffusing in the dissipative non-orthogonal 
environment with various strengths of coupling between the two degrees of freedom is systematically 
calculated. The optimal diffusion path of the particle in a non-Ohmic damping system is revealed to 
have a probability to return to the potential valley under the combined influence of the off-diagonal 
system tensors. 
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I. INTRODUCTION 



The problem of escaping from a metastable state potential is of ubiquitous interest in 
almost all scientific areas in particular the study of nuclear reactions. The related reactive 
system is usually modeled within the framework of one-dimensional (ID) standard Brow- 
nian motion [l-4]. However, since it leaves out the correlation between different degrees 
of freedom while many processes obviously involve more than one degree of freedom, the 
classical ID model fails to describe satisfactorily the dynamical evolving of a real reactive 
process. Basing on these considerations we have recently generalized the ID model to the 
two-dimensional (2D) case by analyzing a set of coupled generalized Langevin eqnarryj^, (J. 

It has been shown that the diffusion in a 2D potential energy surface (PES) includes 
nutritious useful information of the reactive dynamics such as there exists an optimal incident 
angle (or an optimal path) for the particle to obtain its maximum probability to surmount 
the PES barrier. This provides a convenient way to understand many stochastic dynamical 
processes such as the fusion of massive nuclei and even the synthesis of super-heavy elements 
because one can easily estimate the reactive probability of a particle by tracing its footprint 
along the optimal path. 

Moreover, It is also revealed that the non-orthogonality of the PES and the off-diagonal 
system parameters is very important in determining whether a reactive processes can be 
easily accomplished or not. Although this was mentioned in some previous studies 0, Q], 
less effort has been made to give a thorough investigation. In particular, no research has 
based directly on the optimal incident angle as far as we have known. Therefore it is very 
meaningful to seek for more detailed information on this subject. 

In this paper, motivated by the interest of better understanding the 2D reactive dynamics, 
we present a relatively systematic study of the optimal incident angle which enables the 
particle to surmount the barrier with maximum passing probability. Firstly, in Sec. [TT] 
a large number of Langevin calculations are performed in the Ohmic damping case by 
simultaneously varying the off-diagonal term of the system tensors. Secondly, in Sec. IHIl 
influence of the non-orthogonality of the system tensors on the non-Ohmic damping diffusion 
process is discussed where a startling potential valley returning behavior of the optimal path 
is witnessed. Sec. IIVI is a summary of our conclusion in which also the implicit application 
of this work is discussed. 
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II. COMBINED INFLUENCE OF OFF-DIAGONAL SYSTEM TENSORS 



In brief, we begin with the optimal incident angle defined in our previous study by tracing 
the minimum value of critical initial velocity i>° in a 2D Ohmic damping barrier surmounting 
pro cess j^]. It reads 

, _ fm 12 ((3 22 a + UJ22) - m 22 ((3 12 a + m 

0m " arCtan V m 11 F 1 {a)+m Vi F 2 {a) ) W 

where Fj(a)s, (j = 1,2) are the related functions derived in the process of Laplace solving 
the coupled Langevin eqnarray m^Xjit) + (3ijXj(t) + u)ijXj(t) = Expressions of them 

are 

Fi(a) = m 22 a 2 + f3 22 a + lo 22 , 

F 2 (a) = -m 12 a 2 - f3 12 a - u 12 , (2) 

with a the largest analytical root of 

(detm)s 4 + (m n f3 22 + m 22 /3 n - 2m l2 f3 l2 )s 3 
+ (det/3 + niuuj 22 + m 22 Uu - 2m 12 u 12 )s 2 

+ (/3nw 2 2 + ^22^11 - 2/3 12 ui 2 )s + detu = (3) 

the symbols with subscripts such as m 12 are the components of inertia (m^), friction (fy) 
and potential-curvature (cjy) tensors respectively, "det" denotes the determinant of each 
tensor. In the Ohmic damping case, all the system tensors in the Langevin eqnarray can 
be considered as invariable constants and correlations of the two components of the random 
force £i(t) obey the fluctuation-dissipation theorem (£i (£)£/(£')) = ^^' m ife Pkj5{t — t'), where 
ks is the Blotzmann constant and T the temperature. 

Noticing, Eq. ([1]) implicitly contains the off-diagonal term of the system tensors. This 
implies the non-orthogonality of the system tensors is an important factor in determining the 
reactive probability of a particle. A combined control of them on the optimal incident angle 
can be expected. The primary purpose of this work is then to obtain some comprehension 
on this phenomenon. 

In the calculations here and following, we rescale all the variables so that the dimensionless 
unit is used. In Fig. [TJ the optimal incident angle is plotted as a function of various off- 
diagonal parameters (p ofF _ dia ) . From which we can see that the optimal incident angle varies 
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FIG. 1: Optimal incident angle as a function of various off-diagonal parameters (p fF-dia)- When 
one is varying the other two remains zero. Other parameters used here are: m\\ = 1.5, m-22 = 2.0, 
011 = 1-8, /3 22 = 1.2, wii = -2.0, u 22 = 1.5. 




FIG. 2: Optimal incident angle as a function of /?i2 for various ujyi (or myi). Identical system 
parameters are used as those given in Fig. [TJ 
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almost linearly from positive to negative as the increasing of u>u (or P12) when there is no 
influence of other off-diagonal ones, while the effect of m\i is completely on the contrary. 
Supposing the diffusion takes place in a 2D x\ — x<i PES and X\ is the potential valley 
direction, it implies that the non-orthogonality of the potential-curvature (or the friction) 
makes the average diffusion path of the particle to turn toward the negative x<i axis, while 
the asymmetry of the inertia has a opposite affection. This reveals, in a certain reactive 
process, the optimal incident angle may to some extent be controlled by the varying of the 
non-orthogonality of the system tensors. The co-operation of the three off-diagonal system 
tensors may lead to a relatively ideal dynamical reactive path for the particle to surmount 
the barrier. 

In order to get more detailed information, we made a thorough analysis about the influ- 
ence of the non-orthogonality on the reactive dynamics by considering simultaneously the 
varying of all the three off-diagonal components. In Fig. [2j the optimal incident angle is 
plotted as a function of /3i2 for various W12 (or m^) at certain m\2 (or ^12)- A common char- 
acter is found as that the optimal incident angle decreases almost linearly as the increasing 
of Pi2- However the decaying rate is discordance in each subgraph. An important conclusion 
can be made by comparing each curve in Fig. [2] as that for a reactive system with definite 
off-diagonal friction tensor the non-orthogonality of the potential-curvature or symmetry of 
the inertia could help the diffusing particle to obtain a big probability to pass the barrier. 
This is of directive significance in the simulation or experimental operation of many reactive 
processes such as the fusion of massive nuclei because accordingly one could try to select 
symmetric collision and non-orthogonal approximation of the potential to obtain an ultimate 
fusion probability. 

For comparison, in Fig. [31 the optimal incident angle is plotted as a function of uu 
(or m 12 ) for various /? 12 . From which we can see, the optimal incident angle decreases 
(or increases) almost linearly with the increasing of U12 (or m.12) at various rates. But 
noticing in subgraphs Fig. [3] (a) and (c) a large value of (3\2 makes it easy for the optimal 
incident angle to reach zero (the potential valley direction) while in subgraphs Fig. [3] (b) 
and (d) a small /? 12 is appreciated. This reveals, given the potential and inertia asymmetry is 
definite, the influence of the non-orthogonality of friction on the diffusion process is relatively 
complicated. This is comprehensible in the fusion process of massive nuclei. From the view 
point of diffusion induces fusion reaction, the friction of a reactive system relies mostly on the 
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FIG. 3: Optimal incident angle as a function of W12 (or 77112) for various /3i2- Identical system 
parameters are used as those given in Fig. [0 

coupling between the system and the bath environment. In the 2D case we concerned, the 
strength of friction is also restricted by the coupling between two degrees of freedom. Thus 
results in for the diffusing particle a time- dependent or coordinate-dependent environment 
which makes the barrier surmounting process complicate. 



III. POTENTIAL VALLEY RETURNING OF THE OPTIMAL PATH 



In order to have a deep comprehension on this phenomenon, we consider a type of time- 
dependent system friction resulted from the non-Ohmic power spectral density Jij(u) = 
/ o; r ) 5 [9 - 12 ] , where 5 is the power exponent taking values between and 2 in which 
5 = 1.0 corresponds to the Ohmic damping case discussed above. 7^ is the symmetrical 
friction constant tensor, and u r denotes a reference frequency induced in order to ensure the 
components of 7^ to have the dimension of a viscosity at any 5. This kind of non-Ohmic 



damping can describe a large group of anomalous diffusions [13l-ll5|. 



The leading part in the following discussion is still based on the optimal incident angle. In 
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FIG. 4: Non-Ohmic optimal incident angle as a function of the power exponents 5 for various 
effective frictions. The system parameters for each curve are (a) 712 = 0.8, 711 = 1.8, 722 = 1-2; 
(b) 712 = 2.0, 711 = 2.5, 722 = 2.2; (c) 712 = 4.0, 711 = 4.8, 722 = 4.2, with co 12 = -0.5 and 
777,12 = 0.6 except those diagonal components of each tensor as those given in Fig. [TJ 

the non-Ohmic damping case it can be obtained by a simple generalization from Eqs.([I])-(j3J) 
with all the components of friction /Sy = 1,2) replaced by the corresponding Laplacian 
transformation /%[s] = 7yS a-1 [6j. Where 7^ = jijU^ sin' 1 (5tt / '2) is the effective friction 
with each component of 7^ being supposed to be a constant independent of 5. Thus the 
optimal incident angle in the non-Ohmic case reads 



arctan 



m 12 022[s]a + U22) ~ m 2 20i2[s]a + Up) 
m n Fx(a) + m 12 F 2 (a) 



(4) 



implicitly containing the power exponent 5. 

In Fig. HJ the non-Ohmic optimal incident angle $ m is plotted as a function of exponent 8 
at various strengths of effective frictions. In which it is revealed that the non-Ohmic optimal 
incident angle evolves as a non-monotonic function of the exponent 5. Given the effective 
friction is relatively strong (seen line (b) and (c) in Fig. H]for example ) the optimal incident 
angle tends to approach zero as 5 is varying from non-Ohmic region to Ohmic case (5 = 1.0). 
This is an un-trivial behavior because $ m = represents the potential valley direction. So 
an amazing prediction can be made as that the optimal diffusion path will in some case 
return to the potential valley as is always expected in the ID model. It is also distinguished 
from the 2D Ohmic damping case where the optimal diffusion path is usually considered 
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FIG. 5: The maximum stationary passing probability as a function of the power exponents 5 for 
various effective frictions. Identical system parameters are used as those in Fig. [4j 

deviating from the potential valley direction [5j. This can be understood from the view point 
of combined influence of the non-orthogonality of system tensors and the coupling resulted 
2D non-Ohmic damping environment. It is their combination who lead the returning of 2D 
non-Ohmic optimal diffusion path to the potential valley. 

This can also be understood by investigating the stationary barrier passing probability 
of the particle which is usually known as a supplemented complement error function P st = 
^Erfc [— (xi(t))/v / 2cr a ; 1 (t)] of the reactive degree of freedom denoted as X\ here 6j. In fig. 
[5] we plot the maximum value of stationary barrier passing probability P s t, m (it can be got 
by tracing the particle along the optimal path) as a function of the power exponents 5 for 
various effective frictions. From which we can see, at most cases of non-Ohmic damping 
(5 7^ 1.0) the maximum stationary barrier passing probability is smaller than the Ohmic 
damping case (5 = 1.0) except for a narrow range of 1.0 < 5 < 1.4. This is because 
different values of power exponent 5 results in different values of critical initial velocity. 
For example, one can find v*\ s=1 . 8 = 2.6417 > ^| 5=0 . 6 = 2.3477 > v c o \ s =i. = 1.8426, while 
Vq\s=i.i — 1.8251 < 1 5=i.o — 1.8426, calculated by using of the system parameters presented 
in Fig. IH^a) and incident angle $ = 0.258rad. Thus it results in the non-monotonic behavior 
of the maximum stationary barrier passing probability P s t,m and then the potential valley 
returning of the optimal path. 
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Although it is not the primary purpose of this work, determination of the fusion prob- 
ability is very important in the study of fusion dynamics. The potential valley returning 
behavior of the optimal path revealed in present study will provide very useful information 
for the experimental studying of a real reactive process. This is because one can accord- 
ingly select the most appropriated combination of the non-orthogonality of system tensors 
or try to make some possible adjustment on the dissipative environment to obtain a big 
fusion probability. For example, symmetric collision and non-orthogonal approximation of 
the potential in Ohmic damping environment and suitable friction strength in non-Ohmic 
case are appropriated as revealed in our study. 

IV. SUMMARY 

In conclusion of this paper, the barrier passage problem of a particle diffusing in a 2D non- 
orthogonal PES is studied in the framework of statistical Langevin reactive dynamics. In the 
whole range of friction strength from Ohmic to non-Ohmic damping, the optimal incident 
angle of the diffusing particle is found to be greatly influenced by the non-orthogonality 
of the system tensors. A type of potential valley returning behavior of the optimal path 
is witnessed in the 2D non-Ohmic damping environment under the combined influence of 
off-diagonal system tensors and the coupling between the two degrees of freedom. The 
result of this work provides useful information to the study of stochastic dynamical reactive 
processes such as the fusion of massive nuclei and those in connection with the synthesis of 
super-heavy elements. 
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